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ABSTRACT 

Aims. To investigate variability and to model the pulsational behaviour of AGB variables in the intermediate-age LMC cluster NGC 
1846. 

Methods. Our own photometric monitoring has been combined with data from the MACHO archive to detect 22 variables among the 
cluster's AGB stars and to derive pulsation periods. According to the global parameters of the cluster we construct pulsation models 
taking into account the effect of the C/0 ratio on the atmospheric structure. In particular, we have used opacities appropriate for both 
O-rich stars and carbon stars in the pulsation calculations. 

Results. The observed P-L-diagram of NGC 1846 can be fitted using a mass of the AGB stars of about 1.8 Mq. We show that the 
period of pulsation is increased when an AGB star turns into a carbon star. Using the mass on the AGB defined by the pulsational 
behaviour of our sample we derive a cluster age of 1.4 x 10' years. This is the first time the age of a cluster has been derived from the 
variability of its AGB stars. The carbon stars are shown to be a mixture of fundamental and first overtone radial pulsators. 
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1. Introduction 

The cluster NGC 1846 belongs to an intermediate age popula- 
tion in the Large Magellanic Cloud (LMC). In contrast to the 
Milky Way the LMC harbours many clusters with an age around 
1 to 5 Gyr (e.g. Girardi et al. |1995l l. These clusters thus provide 
an interesting possibility to study the late stages of stellar evolu- 
tion for stars around 1.5 to 2.5 Mq, while the globular clusters 
belonging to the Milky Way trace the evolution of stars of only 
up to about 0.9 M0. 

While there is a general agreement in the literature that NGC 
1846 is indeed an intermediate age cluster, the published values 
on global parameters of this system show some scatter Derived 
metallicities can be roughly divided into two groups, one around 
[Fe/H]=-1.5 (Dottori et aL .1983t Leon ardi & Rose i2003i) and 
one around [Fe/H]=-0.7 (Bessell et al.[T983l Bica et al. [T986l 
Olszewski et al. fTWTI Beasley et al. |2002] l. Probably the fir st de- 
termination of the metallicity was done by Cohen ( |1982| l giv- 
ing a value in between the two discussed metallicity levels, 
namely [Fe/H] =- 1.1. The most recent determination of the clus- 
ter's metallicity is given by Grocholski et al. (120061) . who derive 
a value of [Fe/H]= -0.49+0 based on the Ca triplet strength of 
17 individual cluster members measured with the VLT. 

Age values for NGC 1846 scatter between less than 1 Gyr 
(Frantsman 1988^ and 4.3 Gyr (Bica et al.l9Wi. Mackey & 
Broby Nielsen (2007) give cluster ages from 1.5 to 2.5 x 10^ 
years based on a comparison of the colour magnitude diagram 
and theoretical isochrones. The range in age thereby results from 
the usage of two different sets of isochrones. These ages corre- 
spond to masses on the AGB of about 1.3 to 1.8 M©. Mackey 
& Broby Nielsen (2007) also report on the probable existence of 
two populations in NGC 1 846 with similar metallicity but sepa- 
rated in age by about 300 Myr. 



A reddening of Av=0.45mag has been determined by 
Goudfrooij et al. (2006 ). Grocholski et al. (|20061 ) note that NGC 
1846 is probably suffering from differential reddening without 
giving any detailed numbers. Keller & Wood (I2006I I give a some- 
what lower reddening of E(B - y)=0.08, i.e. an Ay value of 
0.25 mag. 

The first studies of the stellar content of NGC 1846 were 
published by Hodge ( 119601 ) and Hesser et al. d 19761 1. A number 
of luminous stars on the Asymptotic Giant Branch (AGB) were 
identified and published with finding charts by Lloyd-Evans 
(1980). Throughout this paper we will use the naming given in 
Lloyd-Evans' publication (LEXX) except for H39, which follows 
the numbering from Hodge ( I1960I ). Details on individual AGB 
stars have been published in a number of papers. Frogel et al. 
( [19801 [T990] l and Aaronson & Mould ( fT985] ) gave near infrared 
photometry and mboi values for most of Lloyd-Evans' AGB stars 
and a few more cluster objects, unfortunately without any find- 
ing chart. Tanabe et al. d 19981 1 searched this cluster for extreme 
infrared stars, i.e. stars with a high circumstellar absorption in 
the visual range, but found none. 

No investigations on the variability of the AGB stars in NGC 
1846 has been published up to now. However, as most AGB 
stars are pulsating (forming the group of long period variables 
or LPVs) it seemed very likely that some light variability would 
be found in these stars. In a previous paper (Lebzelter & Wood 
'20051) we showed that the variability analysis of LPVs in a sin- 
gle stellar population like the globular cluster 47 Tuc allows one 
to investigate various fundamental aspects of AGB stars like the 
evolution of the pulsation mode or mass loss. With NGC 1846 
we have chosen an interesting alternative target since its AGB 
stars are more massive and include a number of carbon stars. 
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2. Observations and data reduction 

For the detection of AGB variables in NGC 1846 and the deter- 
mination of their pulsation period two data sets were used. We 
started a photometric monitoring program using ANDICAM at 
the 1 .3m telescope on CTIO which belongs to the SMARTS con- 
sortiunfl. The plan was to obtain images of the cluster in V and 
/ at a rate of approximately twice a week. However, this plan 
was not feasible due to a strong pressure on the observing time 
at this instrument. Still we got a total of 63 epochs between JD 
2453584 and JD 2454048 with a gap between JD 2453856 and 
JD 2453950. 

We used the pipeline reduced V and / images provided by the 
SMARTS consortium. The detection of variables and the extrac- 
tion of the light curves was done using the image subtraction tool 
ISIS 2.1 (Alard |2000l l. The resulting light curves were searched 
for periodicities using Period98 (Sperl |1998i ). A maximum of 
three periods was used to fit the light change. 

It turned out that the data set obtained at CTIO was very 
valuable to detect periods below 100 days, but that it was too 
short for a reliable determination of the long periods occurring in 
some of the AGB stars of NGC 1846. To improve our period de- 
termination we thus decided to access the MACHO database as 
the cluster was covered by this survey. We used the coordinates 
of the variable AGB stars detected within our CTIO time series 
to extract the corresponding MACHO light curves. Two addi- 
tional long period variables not apparent in the CTIO data were 
also found. The MACHO data consisted of about 250 epochs 
over the time span between JD 2449001 and JD 2451546, i.e. 
about 2500 days. The cluster was in the overlap region of two 
MACHO frames so that for part of the sample two series of mea- 
surements were available. Periods were determined in the same 
way as for the CTIO data. Data both from the blue and from the 
red pass band were analyzed separately and gave very similar 
results. Also, time series from the two different MACHO frames 
resulted in the same periods within a few percent. Due to the sig- 
nificant separation in time and due to the different filter systems 
no attempt was made to combine the MACHO and the CTIO 
data. 

Periods determined from both data sets were then compared 
with each other The agreement for periods shorter than about 
100 days was quite good (differing by a maximum of 10 percent) 
taking into account the significant irregularities found on top of 
the regular light change. An average period was calculated from 
the two data sets and used in the further analysis. For LE16 and a 
previously unidentified variable, LWl, the MACHO light curve 
was not useable. For H39 no counterpart could be identified in 
the MACHO data base. All longer periods were taken from the 
MACHO data. Compatibility of the MACHO period with the 
data set from CTIO was checked by eye in all these cases. 

The uncertainties of the periods, calculated with a least 
square fit approach, are less than a few percent. For the MACHO 
data set, photometric errors were taken from the MACHO 
database. Typical values are below 0^01. For the CTIO data set 
we estimate a typical error of 0^01. 

As part of the CTIO time series we obtained near infrared 
time series in J and K using the same instrument. This data set 
covers roughly three months and allows us to estimate the light 
variability in the near infrared. Standard data reduction was ap- 
plied to the near infrared images using IRAF. The background 
was subtracted using dithering. The J and K magnitudes were 
transformed to the 2MASS system using the transformation from 



' www.astronomy.ohio-state.edu/ANDICAM/ 



3 r LW4 - 

1 r - 



9 H \ 1 1 1 1 1 1 1 1 \ 1 1 1 \ \ \ 1 1 \ \ h 




_J I I I I I I I 1 I I I I , I L_ 

0.6 0.7 0.8 0.9 

X(jum) 



Fig. 1. Spectra of the NGC 1846 AGB variable stars without a 
previous spectral type. 

Carpenter (1200 11 1. A set of reference stars within the cluster with 
near infrared photometry data from 2MASS was chosen to set 
the zero point for each epoch. These stars remained constant 
within +0^02 relative to each other on all frames. This value can 
also be seen as the typical accuracy achieved for the variables. 

Additional J and K photometry was obtained with the near- 
infrared array camera CASPIR (McGregor et al. '1994') on the 
2.3m telescope of the Australian National University at Siding 
Spring Observatory. The J and K magnitudes on the CASPIR 
system were converted to the 2MASS system using the transfor- 
mations in McGregor ( 119941 ) and Carpenter (I2001I I. From these 
data we also obtained near infrared photometry for a large num- 
ber of non variable, low luminosity AGB and RGB stars required 
for calibrating the models (see below). 

For 5 variable stars with a previously unknown spectral type, 
spectra were obtained with the Dual Beam Spectrograph on the 
ANU 2.3m telescope at Siding Spring Observatory. The spec- 
tra are shown in Fig. [T] and the spectral types are listed in Table 
|2] For the 0-rich stars, the spectral types were based on the 
strengths of the TiO bands at 616 nm and 705 nm. LW4 appears 
to be of late K spectral type while LWl, LEI 4 and LEI 5 are of 
early M spectral type. The C2 bandhead at 564 nm and the strong 
CN bandheads at 692 and 788 nm clearly identify LW2 as a car- 
bon star (see Turnshek et al. |1985| for an atlas of K, M, S and C 
star spectra). 

3. Results 

A total of 23 red variables were detected within the 
CTIO/ANDICAM field. Nineteen of th ese co uld be identified 
with AGB stars found by Lloyd-Evans ( I1980I I. For two of the 
stars, LE9 and LEI 8 variability was detected, but no period 
could be determined. Four further variables had no previous 
identifier in the literature. Their coordinates are given in Tab.[T] 
Similar to our approach in Lebzelter & Wood (120051 ). we named 
them with a prefix LW followed by a number. 
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Fig. 2. MACHO light curves of the variables of 
our survey showing long periods. Fourier fits 
with up to two periods are shown. Note that the 
y-axis scale varies from star to star For LEI 3 
only part of the light curve is plotted for a better 
illustration of the short period light change. 



Table 1. Coordinates from the 2MASS catalogue for variables 
with no identifier in the literature. For cross correlation with in- 
frared sources identified by Tanabe et al. (I2004ai) see sect.|5]l. 



Ident 


RA (2000) 


DE (2000) 


comment 


LWl 


05:07:36.96 


-67:27:56.9 




LW2 


05:07:40.26 


-67:27:40.3 




LW3 


05:07:54.46 


-67:24:46.2 


prob. non-member 


LW4 


05:07:34.15 


-67:27:48.5 





For stars with their light change dominated by long peri- 
ods, MACHO light curves with corresponding Fourier fits are 
presented in Fig.|2] For the shorter period stars we show the 
CTIO data with the corresponding Fourier fits in Fig. [3] LW3, 
shown in Fig.[3j is obviously a nicely periodic (83 d), red 
((J - K)2MASS =1.17) variable but, as it is located more than 3.5 
arcmin away from the cluster center, we have to strongly doubt 
its membership to the cluster This star will be excluded from 
the further analysis. Thus we have now a sample of 22 red vari- 
ables associated with NGC 1846, 20 of them showing a periodic 
light change. The determined periods of these stars together with 
other basic data are listed in Tab.|2] 

At a first glance two features of our sample are remark- 
able. First, multiperiodicity seems to be a common phenomenon 
among these stars. Only the weakest stars of our sample could 
be fitted with a single period. Thereby we observe two kinds of 
multiple periods, on the one hand two periods with a ratio of 



about 2, and on the other hand long secondary periods, a well 
known feature of AGB variables in the Galactic field and in the 
MCs (e.g. Wood et al. 2004i Wood l2007b . This long secondary 
period can be the dominant one as can be seen e.g. in LE12 and 
LEO (see Fig.|2]i. Second, there are no large amplitude stars in 
NGC 1846, i.e. the classical miras seem to be missing. The light 
variations in the visual hardly exceed 1 mag. However, we find 
a considerable fraction of stars in the period range typically oc- 
cupied by the miras (i.e. around 300 days). 

Tanabe et al. ( I2004bl l found a rapid increase in light ampli- 
tude once the star reaches the tip of the AGB. We see a similar 
relation between the light amplitude and mean K and J - K val- 
ues (Fig. nil. 

The near infrared measurements show only a slight scatter 
typically of the order of the estimated error of the individual 
measurements. This is not surprising given the low light ampli- 
tudes found in the visual data. One star, LE 5, clearly exhibits 
a light change above the error level. Its light curve is shown in 
Fig H The same behaviour of the light change as in the visual 
can be seen clearly. As the time series spans only a time interval 
of about 3 months long time changes are not properly covered. 
Thus we cannot exclude that the true light amplitude in the near 
infrared is larger than our estimate. The variation of about 0.2 
mag seen for LE5 in J and about 0.15 mag in K correspond to 
about 0.4 mag in the visual light curve. As the light amplitude 
of our sample stars hardly exceeds 1 mag in the visual we can 
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Fig. 3. ANDICAM/CTIO light curves of the 
variables of our survey dominated by short pe- 
riods. Fourier fits with up to two periods are 
shown. Note that the y-axis scale varies from 
star to star. LW3 is probably not a member of 
NGC 1846 (see text). 




Fig. 4. Amplitude vs. K (right panel) and J - K (left panel). No 
correction has been obtained for the difference in the filter pass- 
band between the MACHO and the ANDICAM time series. The 
amplitude given corresponds to the maximum light difference 
found in our time series. O-rich and C-rich stars are marked by 
dots and circled dots, respectively. 



estimate that the total amplitude in J and K will be less than 0.6 
and 0.4 mag, respectively. 

Near infrared colour changes are also very small (see Fig.|5]). 
We used our time series to calculate mean K and J - K values 
for the further analysis. Our values are in reasonable agreement 
with the 2MASS data. 



4. Pulsation models 

Given the large number of LPVs now known in NGC 1846, we 
can now attempt to model the period-luminosity (PL) laws. For 
our models, we adopt the following parameters: distance mod- 
ulus (m - M)i/=18.54 and reddening E{B - y)=0.08 (Keller & 
Wood (2006), helium mass fraction Y=0.27 and a metal abun- 
dance Z=0.004 (following Grocholski et al. |2006] l. 

Below, we derive the mass of the AGB stars by fitting their 
period-luminosity (PL) relations, using the values mentioned 
above as a guide only. We note that the variables in this clus- 
ter are mostly of relatively small amplitude so that high mass 
loss rates are unlikely to have begun. It is also known that red gi- 
ant mass loss prior to the superwind phase in intermediate mass 
stars with M ~ 1 .5 M© is much less significant than in lower 
mass red giant stars such as those in globular clusters. We there- 
fore assume a constant mass along the AGB where our variables 
are found. 

The static and linear non-adiabatic pulsation models were 
created with updated versions of the pulsation codes described 
in Fox & Wood (1982). For temperatures above \ogT - 3.75, 
the interior opacities of Iglesias & Rogers ( 119961 ) were used, in- 
cluding tables with the carbon enhancements required for carbon 
stars. For temperatures below log T - 3.75, a scheme similar to 
that of Marigo ( 2002 ) was used for both oxygen-rich and carbon- 
rich opacities. Molecular contribution due to CN, CO and H2O 
and a contribution due to metals of low ionization potential were 
added to low temperature opacities for zero-metallicity mixtures 
of Alexander & Ferguson (1994). Details will be given else- 
where. The core mass was obtained from the L-M^ relation 
ofWood&Zarro a983i 
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Table 2. The AGB variables of NGC 1846. Reddening corrected 
photometry values are given. Periods are given in days. Multiple 
periods are listed in consecutive lines. A colon marks an un- 
certain period. Spectral types (last column) are from Frogel et 
al. fl990. see references therein) except for LE8 and LEO which 
are taken from Lloyd-Evans ( |1983l l. LW3 is probably not a clus- 
ter member and thus not listed. 



Ident P K J - K log L/Lq log T^g Spec. 

C 

c 
c 
c 
c 

M 

s 

K 
C 

C 

S 

M 

M 

M 
M 
M 
M 
C 
K 



Since we are computing pulsation periods, it is important to 
get the radii (and hence T^s) of the models correct. In order to 
produce the correct T^s on the AGB, the model AGB had to 
be made to coincide with the AGB of NGC 1846. The observa- 
tional AGB was obtained by converting K and J-K of the vari- 
ables and non-variables in the cluster to Mk and {J-K)q using 
the distance modulus and reddening given above. Then, for O- 
rich stars, Mk and iJ-K\) were converted to log L/Lq and log Teff 
using the transforms in Houdashelt et al. ( |2000al |2000bl l (the 
2MASS J and K values were first transformed to the SAAO sys- 
tem, similar to that of Houdashelt et al., using the transforms 
in Carpenter 12001) . For C-rich stars, Mk and {J-K)q were con- 
verted to log L/Lq and log T^ff using the {{J-K)o, T^ff) relation 
from Bessell, Wood & Lloyd Evans (|1983l l and the bolomet- 
ric correction to K from Wood, Bessell & Fox (11983). Figure |6] 
shows the NGC 1846 stars in the theoretical HR-diagram. 

In order to match the models to the observed AGB, the 
mixing-length in the convection theory was set so that the O- 
rich, Z=0.004 models matched the T^g of the O-rich stars in 
NGC 1846. As noted above, all our models have the same mass 
of 1.8 Mq: as explained below, this mass is required in order to 
get the correct pulsation periods. It was found that a constant 
mixing-length gave an AGB that decreased its Tgff too rapidly 
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Fig. 5. Near infrared light variations of LE5 in J (middle panel) 
and K (lower panel) from CTIO/ANDICAM data. For compari- 
son the visual light curve (CTIO/ANDICAM) for the same time 
span is shown (upper panel). 



with increasing L (this is a general problem with all models e.g. 
the tracks or isochrones of Girardi et al. 2000). It was found that 
the ratio of mixing-length to pressure scale height a needed to 
increase slowly with L in order to obtain the correct slope for the 
theoretical giant branch (a was found to vary linearly with L"^). 
This can be seen in Table |3] 

It is clear from Figure |6] that above about logL/Lo = 3.8 
(L = 63IOL0), all stars in NGC 1846 are C stars. For L > 
6OOOL0, we therefore assumed that C/O increased linearly with 
logL/L©: this corresponds to a constant rate of carbon dredge 
up since '""^'^ averaged over a shell flash cycle is constant for 
AGB stars (Wood 1974| . We assumed that C/O reached 3 at 
L - 9000Lq, roughly the luminosity of the brightest AGB star in 
the cluster We have no estimates of C/O for the stars in the clus- 
ter but values of C/O ~ 3 at the end of AGB evolution are indi- 
cated by planetary nebulae in the LMC (Stanghellini et al. 2005). 

The model AGB created with the varying C/O ratio and C- 
rich opacities is shown in Figure |6l As soon as C/O passed 
through 1.0, Teff decreased rapidly as the CN molecule and its 
opacity became more prominent. When C/O ~ 1.3, essentially 
all the N atoms are bound to C atoms which are not bound up 
in CO. Thus the decrease in Teff is not as rapid with further in- 
creases in C/O. 

It is notable that the model Teff values are not as cool as the 
Teff values estimated for the observed stars by converting J-K 
to Teff. As shown below, the model Teff values lead to quite a 
good match between observed and predicted pulsation periods, 
suggesting that the model Teff values are reliable. This would 
indicate that the most luminous and cool C stars in NGC 1846 
are warmer than estimated here. This could be because they have 
significant circumstellar shells leading to reddening of the stars, 
or because the conversion from J-K to Teff is incorrect. 

The dashed line in Figure |6] is the track that would be fol- 
lowed by a C-rich star (with C/O = 2 in the case shown) un- 
dergoing a helium shell flash. During the He-burning part of the 
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Fig. 6. The tip of the NGC 1846 giant branch in the HR-diagram. 
The variable and non-variable stars are shown along with their 
spectral types. The solid line is the theoretical giant branch, 
which has C/O increasing when L is greater than 6000 (see 
text). Values of C/O are indicated at several points on the theo- 
retical track. The dashed line is the line which a C-rich star will 
follow when its luminosity declines and rises during a helium 
shell flash cycle. 

shell flash cycle, AGB stars of mass 1-2 decline in lumi- 
nosity by roughly a factor of 2 ( Vassiliadis & Wood I1993I I and 
spend a reasonable fraction of the cycle time near the low lu- 
minosity part of the cycle. We suggest that the two stars with 
log(L/L0) < 3.7 and log Teff < 3.52 are stars in this phase (both 
stars are C stars). 

5. Discussion 

Having matched the O-rich models to the observed values of 
Teff , we ensure that the models matched the observed pulsation 
periods of the O-rich stars. If not, the mass was altered, and new 
calculations of the giant branch track and the pulsation periods 
were made. It was found that a mass of 1.8 gave the best fit 
between observations and theory for the O-rich stars. This cor- 
responds to a cluster age of 1.4 x lO' years using the isochrones 
of Girardi et al. (20000. The fit is shown in Figure Q and the 
model details are given in Table [3] It is clear that all the O- 
rich stars, with the exception of LE 17 near logP = 1.8 and 
logL/Lo = 3.49, are first or second overtone pulsators (the stars 
with long secondary periods are not considered here). 

It is clear from the models that as the C/O ratio increases 
beyond 1.0, the periods of all modes increase rather rapidly 
due to the increase in radius and decrease in T^s caused by the 
molecules in the C-rich atmospheres. Furthermore, it can be seen 
that the C stars fall quite nicely on two sequences correspond- 
ing to pulsation in the first overtone and fundamental modes. In 

^ Using the BaSTI isochrones (Bedin et al. 12005b we find an age of 
1.9 X 10' years for a 1.8 M© star to reach the AGB (using the same 
parameters as Mackey & Broby Nielsen [2007l >. 
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Fig. 7. The log(L/L0)-log P diagram for NGC 1846 stars and 
models. The fundamental mode and the first three overtones are 
shown for models with L near the H-burning maximum of the He 
shell flash cycle (solid lines). For AGB stars in the He-burning 
phase of the flash cycle, only the fundamental and first overtone 
modes are shown (dashed fines). 

the absence of these C-rich models, it would have been natural 
to explain the C stars as fundamental mode pulsators scattering 
broadly around an extension of the fundamental mode O-rich 
sequence. 

The fit of the models to the observed C star points is not per- 
fect. This is not surprising given the rather elementary method 
used to compute the C-rich opacities, combined with the fact that 
we do not really know the C/O ratios in these stars. Furthermore, 
as mentioned above, the stars move up and down in luminosity 
during a He-shell flash cycle by factors of about 2 in luminosity 
(Vassiliadis & Wood [1993] ). The dashed lines in Figure [T] show 
how the fundamental and first overtone periods would vary over 
a shell flash cycle for a star with C/O = 2. It seems most likely 
that the two C stars with log L/Lq < 3.6 are near the luminosity 
minimum of a shell flash cycle (where AGB stars spend about 
20% of their time). 

In general, the theoretical and observed periods for the C- 
stars agree reasonably well. This indicates that the radii and T^s 
values of the models agree well with the real values for these 
stars and the assumption that the temperatures indicated by the 
near infrared colours are misleading. 

It is somewhat strange that the most luminous of the C stars 
appear to have shorter pulsation periods than the less luminous 
C stars. One way for this to occur would be to decrease the 
C/O ratio in the most luminous C stars by hot-bottom burning. 
However, at 1.8 M©, the mass is too low for hot-bottom burn- 
ing, since AGB masses greater than about 4 M0 are required for 
hot-bottom burning to be significant. Another explanation may 
be that the nonlinear pulsation period differs from the linear pul- 
sation period when the pulsation amplitude increases, as noted 
in Lebzelter & Wood (2005). The visual light amplitudes give 
no indication for this as the amplitudes of the carbon stars are 
all very similar (see Fig.|4]l. In the infrared, as shown in Fig.|5] 
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M/Mg 


Mc/Mg 


m,, 


log Teft 


C/O 


Po 


Pi 


P2 


Ps 


AGB 


stars with L near the H-buming maximum of the flash cycle 








6000 


1.8 


0.596 


1.882 


3.5516 


0.324 


167.7 


88.7 


56.8 


49.1 


6069 


1.8 


0.597 


1.884 


3.5503 


0.400 


171.6 


90.4 


57.9 


50.2 


6256 


1.8 


0.601 


1.890 


3.5490 


0.600 


179.1 


93.5 


59.8 


52.2 


6448 


1.8 


0.604 


1.895 


3.5465 


0.800 


188.8 


97.5 


62.4 


54.8 


6547 


1.8 


0.605 


1.897 


3.5447 


0.900 


195.0 


100.0 


64.1 


56.4 


6647 


1.8 


0.607 


1.900 


3.5409 


1.000 


204.6 


104.0 


66.9 


59.5 


6748 


1.8 


0.609 


1.903 


3.5248 


1.100 


235.4 


117.1 


77.9 


70.0 


6851 


1.8 


0.611 


1.905 


3.5139 


1.200 


260.4 


127.1 


86.8 


76.8 


7062 


1.8 


0.614 


1.911 


3.5025 


1.400 


293.3 


139.4 


98.8 


85.2 


7734 


1.8 


0.626 


1.927 


3.4849 


2.000 


371.4 


166.1 


126.8 


104.0 


9000 


1.8 


0.647 


1.955 


3.4652 


3.000 


523.7 


210.8 


175.7 


140.1 


AGB stars in the He-burning phase 


of the flash cycle 










7734 


1.8 


0.626 


1.927 


3.4849 


2.000 


371.4 


166.1 


126.8 


104.0 


4500 


1.8 


0.626 


1.836 


3.5140 


2.000 


163.3 


89.3 


59.1 


53.4 


3500 


1.8 


0.626 


1.800 


3.5276 


2.000 


114.9 


66.4 


43.6 


31.5 



Notes: All models helium abundance Y = 0.26. The initial metal abundance (before third dredge-up) is Z = 0.004. IfHp is the ratio of 
mixing-length to pressure scale height. Pq, Pi, P2 and P3 are the linear nonadiabatic periods (in days) of the fundamental, 1st, 2nd and 3rd 
overtone modes. 



the star with the largest amphtude in our data set is LE5 which 
is also the most luminous star of our sample. This may indeed 
indicate that the shorter periods at the tip of the AGB can be ex- 
plained by nonlinear effects. However, as mentioned above, our 
near infrared time series are too short for a definite conclusion 
on this point. 

A reduction in stellar mass due to significant amounts of 
mass loss can not explain the shorter periods of the most lumi- 
nous C stars since this would increase the pulsation period. In 
general, we have not found it necessary to invoke any significant 
reduction in mass with luminosity up the AGB in these calcula- 
tions. We note that this contrasts with the stars in 47 Tuc studied 
by Lebzelter & Wood (2005 ) where a significant amount of mass 
loss was found on the AGB. As noted briefly in Section 4, this is 
as expected. The stars being studied here and in 47 Tuc appear to 
have relatively low mass loss rates where a Reimers' mass loss 
law may apply. In such a case, the mass loss rate is proportional 
to LR/M. For a given L, this functional dependence suggests that 
the intermediate mass (~1.8 M©) stars in NGC 1846 will have 
lower mass loss rates than the low mass (~0.9 M0) 47 Tuc stars 
because of the higher mass. In addition, the giant branch temper- 
ature increases with stellar mass so the radius will be smaller at 
a given L in NGC 1846. The consequence of the expected lower 
mass loss rate in NGC 1846 is that a smaller amount of mass is 
lost on the giant branches up to a given L. The fractional change 
in mass in NGC 1 846 is even smaller because of the higher mass. 

An independent check of the mass loss can be obtained from 
mid-IR photometry looking for indications of an infrared excess 
due to dust emission. NGC 1846 was covered by the SAGE mid- 
IR survey (Meixner et al. |2006l ). Blum et al. ( [2U06b find [8]-[24] 
as an indicator for mass loss. For five of the carbon stars in our 
sample (LEI, LE2, LE3, LE6 and LEI 1) the corresponding pho- 
tometry can be found in the SAGE point source data base. All 
of them show an [8]-[24] colour typical for low mass loss stars. 
We recall that Tanabe et al. ( |1998ll2004al) did not find any dust 
enshrouded objects in this cluster. For two sources Tanabe et al. 
(2004a) list no identification from the literature. Their object 16 
is obviously the variable LW2. Object 8, the weakest object in 
their list of mid infrared sources, corresponds to LW4. 



Using the SAGE list of point sources we checked for bright 
objects at 24jum with no or a very weak optical counterpart. 
Indeed we found three sources fulfilling this criterion, but their 
membership to the cluster is not very likely. In the appendix we 
give a more detailed description of these three sources. Thus we 
can derive from our findings that the AGB stars of NGC 1846 
lose their mass in a quite short time at the end of the AGB phase. 

The location of LEI 7 in the P-L-diagram remains somehow 
a mystery. While the two carbon stars found at logL/L© = 3.57 
nicely agree with the expected location of a TP- AGB star close 
to its luminosity minimum, it is unlikely that this explanation 
also holds for LEI 7, which is the O-rich star slightly below in 
the P-L-diagram (Fig.|7j. If it would be in the minimum of its 
TP cycle, it would be in a region exclusively occupied by C-rich 
stars during its maximum. A possibility may be that the star is 
located at the "knee" of the fundamental mode sequence at its 
maximum. Alternatively, the derived temperature and luminos- 
ity of this star may be wrong due to higher reddening of this 
object. Circumstellar reddening probably plays only a minor ef- 
fect as there is no clear indication for this from mid-IR photom- 
etry (SAGE; Tanabe et al. '2004a). We further note that the star 
nicely fits onto the giant branch of NGC 1846 (Fig.|6|), thus there 
is no indication for a significant reddening. For the same reason 
we may safely assume that the star is indeed a member of the 
cluster. Further investigations are required to understand the be- 
haviour of this star. 



6. Summary 

Twenty two long period variable stars have been detected in the 
LMC cluster NGC 1846. Most of them were previously known 
AGB stars in this cluster We have modelled the pulsation of 
O-rich and C-rich AGB stars in NGC 1846, taking into account 
the effects of the C/O ratio on the stellar structure. The stars 
fit nicely on several parallel logP-logL-relations. Similar to 47 
Tuc (Lebzelter & Wood 120051 . the LPVs populate the first and 
second overtone sequence at lower luminosities, while the more 
luminous stars are found on the first overtone and fundamental 
sequence. In addition to the separation by luminosity, there is 
also a difference related to the atmospheric chemistry. The O- 
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rich AGB variables are found to be first and second overtone 
pulsators. We show for the first time that the period of pulsation 
is increased when an AGB star turns into a carbon star because 
of the expansion of the star and the lowering of Tgff. The car- 
bon stars are shown to fall on two sequences in the P-L-diagram, 
corresponding to pulsation in the fundamental and first overtone 
modes. Oxygen rich models would give a very poor fit to the PL 
relation for C stars and would predict fundamental mode pul- 
sation only. We find that there has been no significant mass re- 
duction from mass loss right to the tip of the observed AGB in 
NGC 1846. This is in agreement with mid-IR data. The mass 
of the AGB stars is ~1.8 Mq, corresponding to a cluster age of 
1.4 X 10^ years. This is the first determination of the age of a 
cluster based on the pulsational properties of its AGB variables 
and in excellent agreement with the value derived by Mackey & 
Broby Nielsen ( I2007I I from isochrone fitting. However, the sep- 
aration of the two populations found by the latter authors is too 
small to be detectable in the pulsational behaviour of the AGB 
stars. 
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Appendix A: 3 very red sources in the field of NGC 
1846 

The SAGE images show three sources in the area of NGC 1846 
that are very prominent at 24 /im and have no clear counterpart 
in the visual. Table lA.ll gives coordinates and flux values for the 
three objects. 

All three sources are extended as can be seen on the F814W 
band HST/ACS images available in the HST data archive (pro- 
gram ID: 9891 , PI: Gerard Gilmore) and shown in Figs. lA.ll 
and IA.2I respectively. Source 1 and 2 are located much closer 
to the cluster center than source 3. The location relative to the 
cluster is also visible from Fig. 1 of Mackey & Broby Nielsen 
(120071 where source 3 is located in the right part of the image 
just above the gap and source 2 on the upper left part. None of 
the three sources could be identified with any object in the list of 
Tanabe et al. (1998, 2004a). 

Source 1 has no detectable counterpart at 2.2yum but prob- 
ably at 3.6, 4.5 and 8//m. The energy distribution is plotted in 
Fig. lA.3l The object is of rather circular shape with an approxi- 
mate diameter of 2 arcseconds. 

Source 2 looks elongated and is a little bit difficult to separate 
from a nearby star, probably a RGB star of NGC 1846 according 
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Table A.l. Infrared sources in the field of NGC 1846. 
Coordinates and fluxes are taken from the SAGE catalogue. 




Ident RA(2000) DE (2000) (mjy) 

Source 1 05:07:33.21 -67:27:06.2 1.346 

Source 2 05:07:34.42 -67:27:07.2 0.799 

Source 3 05:07:34.94 -67:29:01.8 0.875 



Fig. A.l. HST/ACS image of source 3 (in the middle above the 
gap). The size of the image is approximately 13x13 arcsecs. 
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Fig.A.3. Mid-IR energy distribution of source 1 from SAGE 
data. 



to its 2MASS J and K values. It is below the SAGE detection 
limit at 3.6 yum, 4.5 fim and 5.8 fim, respectively. A source is de- 
tected at 8 fim at a flux level of 0.585 mJy. The object has a size 
of about 3 to 4 arc seconds. 

Source 3 is of slightly elongated shape with an approximate 
diameter of 5 arc seconds. The image reveals a central concen- 
tration and a few spots possibly related to this source. The central 
source is not present in the 2MASS point source catalogue and 
is also not visible on the 2MASS K frame. This object is most 
likely a background galaxy. Interpretation of sources 1 and 2 is 
more difficult. We suspect that due to their sizes and redness both 
are background galaxies. They seem to be far too weak at wave- 
lengths shortward of 20 fim to represent dust enshrouded objects 
within the cluster. 



